Space and Ground Based Pulsation Data of 77 Bootis Explained 
with Stellar Models Including Turbulence 
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ABSTRACT 

The space telescope MOST is now providing us with extremely accurate low frequency p-mode 
oscillation data for the star ij Boo. We demonstrate in this paper that these data, when combined 
with ground based measurements of the high frequency p-mode spectrum, can be reproduced with 
stellar models that include the effects of turbulence in their outer layers. Without turbulence, 
the / = modes of our models deviate from either the ground based or the space data by about 
1.5-4 /iHz. This discrepancy can be completely removed by including turbulence in the models 
and we can exactly match 12 out of 13 MOST frequencies that we identified as I — modes in 
addition to 13 out of 21 ground based frequencies within their observational 2a tolerances. The 
better agreement between model frequencies and observed ones depends for the most part on the 
turbulent kinetic energy which was taken from a 3D convection simulation for the Sun. 

Subject headings: stars: evolution — stars: individual (r) Boo) — stars: oscillations — turbulence 



1. Introduction 

Until now, high precision stellar modeling of the 
outer stellar layers was needed solely for helioseis- 
mic studies. In order to reproduce the observed so- 
lar p-mode oscillation frequencies, Li et al. (2002) 
demonstrated that the proper inclusion of turbu- 
lence improves the observed solar high frequency 
p- modes from a maximum deviation of 15 /iHz at 
4000 /iHz for a model without turbulence to 5 /iHz 
for a model with turbulence. 

The inclusion of turbulence is a twofold prob- 
lem. It consists of realistically modeling turbu- 
lent convection in the outer layers and then in- 
cluding simulation data in stellar models. Semi- 
analytical models for turbulent convection have 
been proposed by Canuto (1990, 1996). His main 
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idea is to include a full turbulence spectrum. 
Canuto's model has been included in stellar codes 
by Canuto & Mazzitelli (1991) and Canuto, Gold- 
man, & Mazzitelli (1996). The free parameters in 
the semi-analytical model were derived from lab- 
oratory experiments of incompressible convection 
and extrapolated to stellar conditions. Using this 
approach, the superadiabatic peak is much higher 
than that of the standard solar model (SSM), 
while the derived p-modes are closer to the ob- 
served values than those from SSMs (Paterno et al. 
1993). 

In another approach, three-dimensional (3D) 
large eddy simulations of deep compressible con- 
vection have been first performed by Chan & Sofia 
(1989). While it is not possible to resolve the 
full turbulence spectrum in 3D simulations, the 
advantage of these simulations is that they are 
essentially parameter-free, provided that the em- 
ployed subgrid model does not significantly modify 
the properties of the large scale dynamics. In the 
early studies, no account was taken of the radiative 
transfer in the 3D simulation, and therefore the 
stellar models that included a parametrized con- 
vective flux from the simulations showed a larger 



1 



discrepancy with observed solar p-mode frequen- 
cies (Lydon, Fox, & Sofia 1992; Lydon 1993). 

Part of these limitations were overcome by Kim 
et al. (1995, 1996) who employed the diffusion 
approximation for the radiation field but conse- 
quently could not include the optically thin part 
of the SAL. Later, Kim & Chan (1997, 1998) em- 
ployed the Eddington approximation for the radia- 
tion field, included a realistic equation of state and 
radiative opacities. Their simulation spanned 5.5 
pressure scale heights and included all of the SAL. 
Demarque, Guenther, & Kim (1999) mimicked the 
effects of the simulations in calibrated solar mod- 
els by increasing the opacity coefficient k which 
decreased the discrepancy between observed and 
model p-mode frequencies. 

Using different numerical methods to solve the 
convective and radiative equations, Stein & Nord- 
lund (1998) also performed full 3D simulations, 
incorporating LTE radiative transfer and a real- 
istic equation of state. The simulation included 
the entire SAL and spanned a total of 11 pres- 
sure scale heights. Rosenthal et al. (1999) used 
averages of Stein & Nordlund (1998)'s hydrody- 
namical simulations to match the simulation to 
an envelope that was constructed with a standard 
mixing length envelope code. These patched mod- 
els showed a better agreement with the observed 
p-mode frequencies than earlier models. 

In this paper we make use of the recently per- 
formed 3D simulations of fully compressible hydro- 
dynamics by Robinson et al. (2003, 2004). These 
efforts build on the earlier work by Kim & Chan 
(1998). While resolving the SAL and covering 7.4 
pressure scale heights in the vertical domain, they 
yielded results in agreement with Stein & Nord- 
lund (1998). Robinson et al.'s studies showed 
that artifacts of the boundary conditions had af- 
fected the Kim & Chan (1998) 3D simulations and 
they ascertained the resolution and domain sizes 
needed to yield physically realistic results. Us- 
ing averages for the turbulent pressure and turbu- 
lent kinetic energy taken from these simulations 
Li et al. (2002) included these effects on the stel- 
lar structure of the one dimensional models within 
the mixing length theory (MLT) framework. So 
far, the best match to observed solar p-modc fre- 
quencies has been achieved with these methods. 



The space mission MOST 1 (Walker et al. 2003) 
is now providing us with new low frequency p- 
modes for r\ Boo (Guenther et al. 2005). With 
the ground based measurements of the high fre- 
quency part from Kjeldscn et al. (2003), which 
are sensitive to the outer stellar layers of the star, 
we demonstrate that the combined data set can be 
matched within the errors when we include turbu- 
lence in the outer layers of this star. The ground 
based data by Carrier, Eggcnberger, & Bouchy 
(2005) are also discussed. 

2. Stellar Models 
2.1. Turbulence 

For the Sun, Li et al. (2002) have devised a 
method to include the effects of turbulence ob- 
tained from 3D hydrodynamic simulations (Robin- 
son et al. 2003) by including both the turbulent 
pressure and the turbulent kinetic energy into the 
ID stellar model within the framework of standard 
MLT. This method produces p-mode frequencies 
that match the measured solar spectrum better 
than an SSM without turbulence. 

We slightly modify the techniques of Li et al. to 
enable us to apply the 3D turbulence data for the 
Sun to our model of rj Boo. Until we have com- 
pleted a full 3D simulation for the outer layers of 
rj Boo, we make use of the solar data and appro- 
priately shift the data to apply it at the correct 
depth in r\ Boo. This shifting is motivated by an 
expected characteristic found in all our 3D simula- 
tions: namely that the SAL peak closely coincides 
with the turbulent pressure peak. 

This is verified for 3D simulations of four dif- 
ferent evolutionary stages of the Sun (ZAMS, cur- 
rent Sun, subgiant and giant). In each case the 
peak of the turbulent pressure closely coincides 
with the peak of the SAL. Figure 1 illustrates 
this property for the Sun. The offset between the 
peaks measured in pressure difference is always 
smaller than Alog 10 (P [dyn/cm 2 ]) =0.1. There- 
fore, in order to apply the solar turbulence data 
to Tj Boo we match the turbulent pressure peak 
from the 3D simulation data with the SAL of 



1 MOST (Microvariability & Oscillations of STars) is a Cana- 
dian Space Agency mission, jointly operated by Dyna- 
con Inc., the University of Toronto Institute for Aerospace 
Studies and the University of British Columbia, with the 
assistance of the University of Vienna. 
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Fig. 1. — Averaged quantities taken from a 3D 
turbulent convection simulation for the Sun. The 
peak of the SAL coincides with the peak of turbu- 
lent pressure. The inset magnifies the peak loca- 
tions to show their relative positions. 

the one dimensional model with a small offset of 
Alog 10 (P [dyn/cm 2 ]) = 0.0758, a mean value de- 
rived from the four solar models. This ensures that 
the solar turbulence data is applied at the correct 
depth of the r/ Boo model. 

T) Boo exhibits a slightly higher effective tem- 
perature of T c ff = 6028 K which amounts to a 18% 
higher flux compared to the Sun. The surface 
gravity of r) Boo is a factor of four lower. Both 
differences can change the relative strength of the 
turbulent pressure and turbulent kinetic energy in 
the 3D simulation but we do not have enough sim- 
ulations available yet in order to extract sensible 
scaling relations. Therefore, no further scaling has 
been performed. 

The refined treatment of the outer stellar layers 
has been implemented in the Yale Stellar Evolu- 
tion Code (YREC). The numerical methods and 
main physics included are described by Pinson- 
neault (1988) and Guenther et al. (1992). The 
most recent improvements other than the inclu- 
sion of turbulence in the outer layers include new 
updates to the equation of state (OPAL 2001 EOS, 
Rogers 2001). 

The high spatial resolution needed for the solar 
models is also required for the r\ Boo models. A 
model typically consists of about 4500 grid points 
which are distributed in order to give smooth runs 
of all variables. In contrast to the solar models we 



find that a more stringent time-stepping is needed 
for r\ Boo in order to yield resolution independent 
p-mode frequencies. This is most likely due to 
•q Boo's more advanced evolutionary stage. This 
demands us to advance the model through at least 
2500 time steps. 

2.2. Starting Model 

Our model construction starts from a model 
for r\ Boo that has been selected with the quanti- 
fied dense grid method (QDG) developed by Guen- 
ther & Brown (2004). The search performed for 
r] Boo is described in detail by Guenther et al. 
(2005). Their best fit model is selected from 
an extended search with different input param- 
eters for hydrogen X = (0.69,0.71), metallic- 
ity Z = (0.02,0.03,0.04) and stellar masses be- 
tween 1.4 M© and 1.9 Mq with a fine grid reso- 
lution of 0.005 Mq. Along each evolution track 
of the models, p-mode frequency spectra for the 
I = 0, 1, 2, 3 modes have been calculated from ra- 
dial order n — 1 to the acoustic cutoff frequency. 

A total of 3 x 10 7 model frequency spectra 
have been compared with eight selected MOST 
frequency measurements that were judged to be 
the most likely members of the I — p-mode se- 
quence. The agreement between model spectra 
and observation is ascertained with the x 2 formu- 
lation (Guenther & Brown 2004). The best model 
consists of a mass of 1.71 ± 0.05 M , (X,Z) = 
(0.71, 0.04), a mixing length of 1.8 and no element 
diffusion 2 at an evolution age of 2.40 ± 0.03 Gyrs. 
This best model with a \ 2 < 1-4 was constrained 
only by the 8 MOST p-modes. No other con- 
straints, such as composition, surface tempera- 
ture, or luminosity were used. Regardless, their 
model which best fits the oscillation data also lies 
within 1(7 of the observationally derived effective 
temperature, luminosity, and metal abundance. 

A new interferometric measurement of 77 Boo's 
radius is now available (Thevenin et al. 2005), 
yielding a radius of R/R Q = 2.68 ± 0.05. The 
best model selected with the QDG search tech- 
nique is fully consistent with this value, since it 
possesses a radius of R/Rq — 2.6842. Guenther 
et al. (2005) show that there is no other model 
within the searched parameter-space that fits both 

2 See Guenther (2004) for an discussion about diffusion in 
rj Boo. 
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the MOST data and, in addition, the effective tem- 
perature and luminosity as derived observationally 
by Di Mauro et al. (2003). Thus, the new inter- 
ferometric radius measurement does not give any 
additional constraint to our modeling, neverthe- 
less it is an essential requirement that our mod- 
els are consistent with this observationally deter- 
mined radius. 

The pulsation spectra computed for this paper 
are calculated from a model with exactly the same 
input parameters but using a slightly different ver- 
sion of YREC that incorporates the newer OPAL 
2001 EOS. For an evolutionary age of 2.409 Gyrs 
we achieve a favorable fit of x 2 < 1-0. This model 
is shown in an echelle diagram in Figure 2 (green 
triangles). All x 2 -numbers are calculated with an 
adopted model uncertainty of 0.05 /iHz, and with 
the exact la uncertainty as quoted by the authors 
for their individual measurements. 

It is important to note that none of the models 
without turbulence - searched for giving a good 
X 2 -fit for eight selected MOST low p-mode fre- 
quencies - match the I = ground based frequency 
measurements satisfactorily. The model frequen- 
cies appear to be at slightly higher folded frequen- 
cies with a difference that increases from 1.5 /iHz 
at 600 /iHz to 4 /iHz at 900 /iHz (Figure 2). These 
differences between model and observations are 
marginally within observational uncertainty up to 
700 /iHz, but above 700 /iHz they are significant. 

2.3. Model Calibration 



1000 -» 



800 



600 



sr 400 



200 - 



1000 



800 - 



600 - 



sr 400 - 



200 - 



1000 - 



800 



600 - 



sr 400 



200 (* 



l=0 



\m 

hah hh ■ 

hAH« * 

hA- ^h& A 

hSBH a 

■ A 

• a hAH 
■ AJj»H 
WH 

HCH 
h«H 

mi 



Measurements 
hen MOST 
H#H MOST (QDG) 
h-BH Kjeldsen 2003 
hAH Carrier 2005 

Models 
■ turbulence 
A no turbulence 



I I I I 



I I I I 



I I I I 



t— r- 



1=1 



hAH hH^H 
hA3H ■ l-EH 
A 

WW h&H 
hAH I-B-* A 

hAH ■ A 
hAH "A 
HSH □ H— I 



*©H hW 

he* 

i — r> 1 1 / rw. r -i 
■LJ I r^^J^ 



I . 1 I I 



I I I I 



I I I I 



■|=2 F 

1 h-AH ■ 

hAH t-am 

— hAH HEU 

I ^ | A 
hAH *&HA 



i h«H 

I h©H A 

• hen 



10 20 30 

Frequency modulo 40.0 [iiHz] 



40 



Next we construct a neighboring model to our 
previously defined best fit model that yields the 
same good agreement to the observed MOST p- 
mode frequencies and in addition includes the ef- 
fects of turbulence in the outer layers. There is 
a fundamental difference in calibrating models to 
rj Boo compared to calibrating models to the Sun. 
For the latter we know the age to high precision. 
Therefore, a solar model is calibrated by evolv- 
ing the model to the exact same age and chang- 
ing two unknown stellar parameters, i.e., mixing 
length and hydrogen mass fraction. To first order, 
the luminosity depends on the hydrogen mass frac- 
tion and the effective temperature is most sensi- 
tive to the mixing length parameter. By attempt- 
ing to follow a similar procedure with rj Boo we 
are faced with the difficulty that neither the age 
nor the mass of this star are known. In order to 



Fig. 2. — Echelle diagram showing the non- 
adiabatic p-mode frequencies derived from a best- 
fit theoretical model without turbulence (trian- 
gles) in comparison to a model with turbulence 
(squares) on top of the ground and space based 
observational measurements. The eight selected 
MOST measurements used in the QDG search are 
marked with filled circles. The quoted error bars 
correspond to a 2 a deviation. 



fit a specific locus in the HRD we can, e.g., hold 
the mixing length constant and only vary the hy- 
drogen content and the age. However, this choice 
is arbitrary and we could with equal justification 
have held either the hydrogen content or the age 
constant, while varying the remaining other two. 
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To find a proper calibration method we look at 
the effect on the p-mode frequencies of changing 
one of the three free parameters (mixing length, 
hydrogen content or age) while keeping the other 
two fixed. Changing the age simultaneously alters 
all frequencies, thus altering the age would destroy 
our good match with the MOST data. Turning 
this finding around, we can view the QDG search 
for a match of the lower frequency p-modes as a 
method for finding the age and locus in the HRD 
of t) Boo. This finding is supported by the more 
rigorous analysis made in Guenther et al. (2005) 
where it is shown that the low frequency p-modes 
anchor the interior structure, hence mass and age, 
effectively. 

With the age and mass being fixed by the low 
frequency p-modes, we conclude that the calibra- 
tion of a model with turbulence has to be per- 
formed the same way as for the Sun by chang- 
ing the mixing length parameter and the hydro- 
gen mass fraction. However, the age parameter 
could be used to fine-tune and improve the x 2 -fit 
of the combined low and high- frequency data sets, 
a possibility not taken advantage of here. 

3. Results 

We now put together the different elements dis- 
cussed in the previous sections in order to derive 
the pulsation spectrum of a model for r\ Boo that 
includes turbulence. The mass, metallicity, age, 
luminosity and effective temperature of r\ Boo is 
derived with the QDG search technique outlined in 
Section 2.2 yielding a model in the subgiant evolu- 
tionary phase at age 2.409 Gyrs, mass of 1.710 Mq, 
metallicity of Z = 0.04 and mixing length of 
1.8 that fits eight selected p-mode observations 
of MOST with x 2 < 1-0. The non-adiabatic p- 
mode frequencies for this model are calculated 
with Guenther's pulsation code JIG (Guenther 
1994). 

Six linear non-adiabatic equations are solved 
which take only into account radiative losses and 
gains. The convective flux is "frozen" out of the 
pulsation equations (see Pesnell 1990, for a de- 
scription on the various ways in which this can 
be done) thus the coupling of convection and the 
oscillations is not accounted for. The calculated 
frequencies are shown in Figure 2. Only one of 
the eight frequencies in the high frequency regime 



matches the ground based data points reported by 
Kjeldsen et al. (2003), on average the model yields 
folded frequencies about 3 /xHz larger. 

Next we include the effects of turbulence in our 
model according to Section 2.1. The free param- 
eters mixing length and hydrogen abundance are 
adjusted slightly to ensure that the luminosity and 
effective temperature of the model with turbulence 
matched the luminosity and effective temperature 
of the model without turbulence. The calibration 
procedure is stopped after luminosity and effec- 
tive temperature match with a relative difference 
better than 5 x 10~ 5 . 

3.1. Radial modes 

Finally, the non-adiabatic I = 0,1,2 p-mode 
frequency spectrum is calculated with JIG and we 
plot the results in an echelle diagram (Figure 2). 
We can see in this figure that the model with 
turbulence still matches the I = low frequency 
MOST data points as required by our calibration 
technique while in addition it reproduces six out of 
the eight ground based I = frequency data points 
by Kjeldsen, five within their la uncertainty and 
one within 2a. Also note that one frequency not 
matching the data is still a match within 3<r and 
that we are using the errors as quoted by Kjeldsen 
ct al. (2003). 

The region of 600 — 650 /iHz where the mod- 
els coincide with the two MOST modes that 
have been independently confirmed by the ground 
based measured modes of both Kjeldsen et al. 
(2003) and Carrier et al. (2005) add credibility 
to our modeling. As already noted in Guenther 
et al. (2005), MOST had measured two modes 
below 200 ^Hz that also fit into the I = sequence 
of our models. 

To provide a more quantitative measure we cal- 
culate the x 2 numbers for the combined data sets 
(MOST plus Kjeldsen et al.) of all I = modes in 
the range 200 — 900 ^Hz. For the model without 
turbulence we get x 2 = 18 compared to x 2 = 2.5 
for the model that includes turbulence. 

In Guenther et al. (2005) the best fit to the 
combined MOST and Kjeldsen et al. modes, again 
only constrained by the oscillation frequencies, 
yielded a x 2 = 2.3. But, importantly, the model 
corresponding to this fit to the oscillation data, did 
not fit rj Boo's location in the HRD. By including 
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turbulence in our model, we fit the MOST oscil- 
lation data, the Kjeldsen et al. oscillation data, 
and the observed position in the HRD (Di Mauro 
et al. 2003; Thevenin et al. 2005) simultaneously. 

Our model with turbulence fits the Kjeldsen 
data much better than the data from Carrier. The 
combined set of MOST plus Carrier data gives a 
X 2 = 18. Since the Carrier data appears consis- 
tently at lower folded frequencies, the standard 
model without turbulence is very far off with a 
X 2 = 131. Hence the model with turbulence is 
still much closer to the Carrier data than a model 
without turbulence. 

The structural difference in sound speed be- 
tween the model including turbulence and with- 
out turbulence is shown in Figure 3 (top panel). 
As expected, the largest deviation of ~ 4% is seen 
within the peak of the SAL (bottom panel). Also, 
the deeper convective layers are affected by about 
1%. The structural differences vanish at pressures 
greater than 10 11 dyn/cm 2 , where the layers are 
fully radiative. 

3.2. Nonradial modes 

Additional information is present in the I = 1 , 2 
modes. From our models we expect some regular 
spacings for the high frequency p-modes, but due 
to t] Boo's advanced evolutionary stage no regular 
spacing should be seen for lower frequencies which 
are subject to strong mode bumping. Unfortu- 
nately, MOST did not see many high frequency 
modes and we must therefore mostly rely on the 
ground based data set. 

We can see from Figure 2 that none of the five 
I = 2 modes from Kjeldsen are matched within 
2a with the model lacking turbulence whereas 
four matches are achieved with the model that 
includes turbulence. The only mode that can- 
not be matched to the Kjeldsen data is one with 
frequency greater than 950 /iHz. Interestingly, we 
cannot reproduce any of the three measurements 
above this threshold (1=0,1,2) suggesting there is 
still room for improvement in our models. We 
identify one additional MOST mode at 648 /iHz 
that fits smoothly into the high frequency I = 2 
sequence. This mode is also better matched by 
the model including turbulence. 

The matches of the model with turbulence to 
the I = 1 ground based data are not as good as for 
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Fig. 3. — Relative difference in sound speed of 
the outer layers of r\ Boo between a model in- 
cluding and lacking the effects of turbulence (top 
panel). The SAL of the model including turbu- 
lence is shown for comparison (bottom panel). 

1 = 2. The model including turbulence matches 
three out of eight measured frequencies in com- 
parison to only two by the model without turbu- 
lence. This ratio is slightly enhanced when we 
include one additional mode found with MOST 
at 828 ^Hz. The more successful fit of the I = 2 
modes can be explained with the finding of Guen- 
ther et al. (2005) who demonstrated that a slight 
perturbation in the model mass of 0.005 M© leads 
only to minimal changes in the / = modes but 
to extreme changes in the I = 1,2 modes with 
particularly large impact on the I = 1 modes. 

Of the many modes smaller than 500 ^Hz that 
are present in the MOST data we plot only those 
which lie near modes predicted by the models. 
Although there are many modes present in the 
MOST data very few modes are actually matched 
by the models and no statistical advantage of 
the models including turbulence can be inferred. 
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Since observations show no evidence for any evenly 
spaced sequences in this frequency domain and 
models predict mode bumping to occur as a re- 
sult of mixed (/-modes for frequencies smaller than 
350 yuHz and mixed p-modes within 350 — 600 pHz 
we are simply not in the position yet to make use 
of this information in our models. 

Of the nine nonradial modes measured by Car- 
rier, none is matched by the standard model fre- 
quencies and only one is matched by a mode from 
the model that includes turbulence. As noted be- 
fore for the radial data, the non radial I = 1,2 
appears at lower folded frequencies. Again, mod- 
els with turbulence fit better to the Carrier data 
than models without turbulence but in cither case 
there remains a large discrepancy between these 
observations and the models. 

Finally, we list the frequencies from our the- 
oretical models in Table 1 together with the ob- 
served frequencies. The observed modes are iden- 
tified with radial orders that most closely match 
our theoretical models. Most of the measurements 
from the ground are identified with a radial order 
n one higher than those in this work. The mixed 
mode character of all modes is indicated by the 
number n g . 

3.3. Origin of improved p-mode frequen- 
cies 

The improvements of the fit between observed 
and model high frequency p-modes arise from 
structural changes in the superficial layers of the 
star, mainly in the SAL of the convection zone. 
However, it is important to note that the form of 
the structure change that reproduces the correct 
shift of the high frequency p-modes is not unique. 
For the Sun, Monteiro, Christensen-Dalsgaard, & 
Thompson (1996) demonstrated that models with 
a steeper and narrower SAL compared to standard 
MLT models can produce a frequency shift that 
brings these models in accord with the observed p- 
mode frequencies. In contrast to this, 3D simula- 
tions for the Sun (Nordlund & Stein 1999) produce 
an SAL stratification that is very close to standard 
MLT. The frequencies of high order p-modes are 
also predicted smaller when calculated from the 
average structure of the 3D simulations, hence, a 
better fit to observations is achieved. Nordlund & 
Stein (1999) attribute this frequency shift in their 
models to the turbulent pressure support and, in 



addition, to 3D effects arising from the net effect 
of the fluctuations of the opacity. 

In the following, we explore how the inclusion 
of turbulence as performed in this paper changes 
the surface layer structure of n Boo and we try to 
identify the characteristic features that lead to the 
correct p-mode frequency shifts. A similar analy- 
sis has been given by Li et al. (2002) for the Sun, 
here, we extend this analysis to n Boo and add 
some more information about the role played by 
the shape of the SAL. 

3.3.1. Turbulent Pressure 

As described in Section 2.1, we account for the 
effects of turbulence on the stellar structure by 
including the turbulent pressure and turbulent ki- 
netic energy taken from a 3D simulation for the 
Sun. Thus, we are able to explore the relative im- 
portance of both effects to the correct shifting of 
the high frequency p-modes. In order to do so, we 
calculate one model for 77 Boo where we include the 
turbulent pressure alone and completely omit the 
turbulent kinetic energy This model is calibrated 
to give the same luminosity and effective temper- 
ature and hence radius as our previous models. 
In Figure 4 we compare this model with turbulent 
pressure alone to the standard MLT model and the 
model including the effects of both turbulent pres- 
sure and turbulent kinetic energy. As can be seen 
in Figure 4, the model with turbulence alone shifts 
the frequencies to higher folded frequencies, hence 
it increases the discrepancy between the model 
and the observations. This finding is supported 
by Li et al. (2002, Fig. 11), where the same effect 
is seen in the case of the Sun. The same behavior 
has been found for the Sun by Balmforth (1992, 
Table 1). 

3.3.2. Turbulent Kinetic Energy 

It is obvious from our model with turbulence 
alone, which fails to shift the p-mode frequencies 
into the right direction, that the main ingredient 
for a better match with observations is achieved 
by the effects of the turbulent kinetic energy To 
show this even more clearly we have calculated 
one additional model in which we artificially in- 
creased the turbulent kinetic energy by a factor of 
two. Again, this model was properly calibrated. 
As can bee seen in Figure 5, increasing the turbu- 
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Fig. 4. — Echelle diagram for a calibrated 
i] Boo model where only the effects of turbulent 
pressure are included (circles) in comparison to 
the standard MLT model (triangles) and our best 
model including the effects of turbulent pressure 
and turbulent kinetic energy (squares). 

lent kinetic energy shifts the frequencies to lower 
folded frequencies in the echelle diagram. 

The plot also demonstrates that a quantitative 
match between observations and models depends 
crucially on the exact magnitude of the turbu- 
lent kinetic energy which we can only derive from 
a complete 3D simulation of the outer layers of 
rj Boo . Since we achieve a good fit to the p- 
mode observations by applying to our ID model 
the effects of turbulent kinetic energy as derived 
from a 3D simulation for the Sun, it remains to be 
shown that a 3D simulation of the outer layers of 
r\ Boo yields comparable values for the turbulent 
kinetic energy. 



Fig. 5. — Echelle diagram for a calibrated 
i] Boo model where the turbulent kinetic energy 
is doubled (circles) in comparison to the stan- 
dard MLT model (triangles) and our best model 
(squares). 

3.3.3. Shape of SAL 

If the effects of turbulent kinetic energy can 
shift the model p-mode frequencies towards the 
observed frequencies, it is illuminating to find out 
about the structural changes that are induced by 
the turbulent kinetic energy onto the outer layers 
of rj Boo. These changes are best seen in the SAL 
(Figure 6). The inclusion of turbulence has the 
effect of shifting the peak of the SAL into deeper 
layers of the stellar envelope. This is the main 
effect responsible for the frequency shift. Also, 
the superadiabaticity is increased (increased peak 
height), but this effect is small. It is worth notic- 
ing that the shape of the SAL is preserved, with 
same half- maximum in all models. In order to 
demonstrate this, we have shifted the SAL of two 
models artificially to make their peak location co- 
incide with the standard MLT model (Figure 7). 
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Fig. 6. — The SAL for three different models: 
a standard MLT model (solid line), our fiducial 
model with turbulent pressure and turbulent ki- 
netic energy included (dashed line) and a model 
with artificially increased turbulent kinetic energy 
(dotted line). 
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Fig. 7. — Same plot as Figure 6 but the models in- 
cluding turbulent kinetic energy have been shifted 
to coincide with the peak of the standard MLT 
model to allow for comparison of the SAL shape. 

We conclude that the inclusion of turbulent kinetic 
energy shifts the SAL peaks into deeper layers of 
the stellar envelope while the shape of the SAL 
is preserved. The deeper SAL location changes 
the run of sound speed in the outer layers lead- 
ing to high order p-mode frequency shifts towards 
lower folded frequencies as required by observa- 
tional data for the Sun and also for r\ Boo. For 
the latter, the exact same amount of turbulent ki- 



netic energy as derived from a full 3D simulation 
for the Sun is able to bring the models in accord 
with the high frequency p-mode observations. 

4. Discussion and Conclusions 

This paper demonstrates that the measured p- 
mode frequencies of n Boo from space (by MOST) 
and the ground (by Kjeldsen) can be jointly 
matched with our theoretical models by includ- 
ing the effects of turbulence in the outer stellar 
layers. We are able to report on a better match 
between theory and observation for a star other 
than the Sun when the outer stellar layers are 
corrected by the effects of turbulence. It is an im- 
portant assumption of this paper, that the effects 
of turbulence on the outer layer of rj Boo can be 
extracted from a 3D hydrodynamical simulation 
of the surface layers of the Sun. 

The high order p-mode frequency shift, that 
brings our model in better agreement with the 
observations, is shown to be a direct consequence 
of the inclusion of turbulent kinetic energy. Al- 
though the turbulent kinetic energy must be 
present in the outer convection zone, it is usually 
disregarded in traditional stellar modeling. The 
turbulent kinetic energy in this study was taken 
from a 3D hydrodynamical convection simulations 
of the Sun, therefore, only a 3D simulation for 
i] Boo can add the final proof that the amount of 
turbulent kinetic energy we use in this paper is 
correct. 

Specifically, we show that the inclusion of tur- 
bulence into our ID stellar evolution models, as 
derived from the solar 3D simulation data, can 
account for the difference of 1.5 — 4/LtHz in the 
echelle diagram between the ground based I = 
data points and the models without turbulence. 
The quantitative agreement for the model includ- 
ing turbulence is excellent for the combined data 
set (MOST plus Kjeldsen) of the I = p-mode fre- 
quencies between 200 — 900 /xHz with a x 2 = 2.5. 

The better agreement between the observa- 
tions and the models that include turbulence, is 
strengthened by comparing the observed ground 
based I = 1, 2 modes. The model with turbulence 
matches four out of nine I = 1 and five out of six 
I = 2 modes. If we combine all 22 observed ground 
based modes with 14 selected MOST modes, a 
model with turbulence reproduces 27 modes versus 
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14 modes for a model without turbulence (within 
2cr error bars). 

MOST sees more spectrum peaks than shown 
in our echelle diagrams. The models predict, that 
we should be seeing mixed modes. However, the 
stellar origin has to be established with more cer- 
tainty. This will be possible with the scheduled 
re-observation of 77 Boo by MOST in 2005. If some 
of the recurring peaks can be identified as nonra- 
dial modes, this will give us excellent additional 
information to test our models. 

Because our theoretical models match the 
Kjeldsen data, they are not a good match to the 
Carrier data, since both data sets do not over- 
lap for the majority of modes which is true for 
both radial and nonradial modes. Regardless of 
this discrepancy, frequencies from models that in- 
clude turbulence come much closer to the Carrier 
measurements. 

Certainly the most important advancement for 
this study is the inclusion of the turbulence data 
derived from a full 3D turbulence simulation for 
77 Boo. We are currently undertaking this task 
and will report on it in the future. Also, the QDG 
search - already being very extensive - must be ex- 
panded in the searched parameter space of hydro- 
gen content, metallicity and convective core over- 
shoot. 

The refinements of stellar evolution theory with 
regard to turbulent convection in the outer layers 
has previously been motivated by attempts to gain 
a better match between theory and observation for 
the frequency spectrum of the Sun. This study in- 
dicates that these refinements to the theoretical 
models are also important for interpreting the ob- 
servational data of r] Boo. We believe that the 
need for extremely precise theoretical models will 
continue to grow as more and more observational 
measurements become available within the young 
field of asteroseismology. 
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Table 1 

Model Frequencies and Observed Frequencies 



Model Frequencies Observations 
Order n p n g With Turbulence Without Turbulence MOST 2005 a ' 6 Kjeldsen 2003 Carrier 2005° 



1 





127.17 


127.35 


126.66/127.91 


2 





170.08 


170.31 


171. 32 d 


3 


(} 


210.96 


211.18 


210. 56 d 


4 





251 .95 


252.19 


nri "70^ 

zoi . / y 


5 





291.78 


292.04 


292. 25 d 


6 





332.45 


332.73 


333. 17 d 


7 





373.02 


373.40 


373. 20 d 


8 





413.60 


414.14 


414. 01 d 


9 





452.51 


453.40 


453.13 d 


10 





490.78 


491.96 


492. 92 d 


11 





529.73 


531.11 


533.0° 


12 





570.29 


571.84 




13 





610.78 


612.71 


610.55 611.0 ±0.5° 610.6° 


14 





650.87 


653.25 


650.37 651.2 ±0.6° 


15 





690.51 


693.21 


690.8 ±0.6° 691.3° 


16 





730.78 


733.54 


732.6 ± 0.4° 729.5° 


17 





770.76 


774.30 


769.4° 


18 





813.18 


816.05 


813.1 ±0.4° 809.2° 


19 





854.09 


857.15 


853.6 ± 0.3° 


20 





894.48 


897.71 


894.2 ±0.6° 891.6° 


21 





933.80 


937.09 




22 





971.60 


1000.99 


974.5 ±0.7° 971.9° 


1 


20/22 


116.00 


104.19 




1 


18 




125.92 




2 


17 


137.60 






2 


16 


142.30 


140.52 




2 


12/10 


177.76 


177.22 




3 


12 




185.19 




3 


11 




196.86 




3 


10 


214.53 






3 


9 


221.96 


221.40 




4 


6/8 


261.27 


259.47 




5 


7 


298.48 


294.38 




5 


6 


304.40 


303.45 




6 


4/6 


339.66 


336.48 




6 


5 




345.58 




7 


3 




379.58 




8 


2 


380.22 






8 


5 


404.03 






8 


4/2 


422.43 


420.77 




9 


2 


448.69 


447.84 




10 


2 


477.37 


476.78 




10 


3 




503.67 




11 


4 


507.43 




512.2 


11 


3 


524.14 


521.11 




12 


1 


553.78 


554.48 


550.3 


13 


1/3 


592.10 


593.18 


589.9 


14 


3 


630.64 


631.91 


629.4 ± 0.3 625.7 


14 


2 


664.90 


658.58 


670.1 ±0.5° 665.4/669. 


15 


2 


678.77 


676.75 




16 


2 


711.75 


713.97 


711.8 ± 0.4 


17 


2 


751.13 


753.68 


749.3/753.4 ± 0.5 748.5 


18 


2 




795.64 


793.1 ±0.7 787.4 


19 


2 


832.52 


835.22 




20 


2 


872.54 


875.18 




20 


1 


911.33 


913.13 
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Table 1 — Continued 









Model Frequencies 


Observations 




Order 


rip 


n g 


With Turbulence Without Turbulence 


MOST 2005 a ' 6 Kjcldsen 2003 


Carrier 2005° 



21 


1 


940.11 


927.53 






22 


1 


955.52 


956.82 


955.6 ± 0.8 


947.6 


23 


1 


994.58 


997.34 






24 


1 


1030.27 


1033.29 


1034.3 ± 0.7 




1-2 1 


30 




129.65 






2 


25 




157.49 






2 


23 




219.36 






3 


19 




200.90 






3 


18 


208.41 


206.27 






3 


16 




224.96 






4 


16 


237.23 








4 


15/14 


246.81 


250.28 






5 


13 


285.04 


283.02 






5 


12 


290.29 


286.99 






5 


11 


310.59 








6 


11 




322.63 






6 


10 


331.04 


328.29 






6 


9 


354.29 








7 


9 


368.25 


368.51 






8 


8 


409.19 


409.66 






9 


7 


448.62 


449.46 






9 


6 




477.88 






10 


7 


486.02 








10 


6 


487.51 


488.28 






11 


6 


526.09 


527.45 






11 


5 


554.59 








13 


5 


607.32 


609.20 


608.1 ± 0.4 




14 


4 


647.83 


649.90 






15 


4 


687.22 


689.95 






16 


4 


727.43 


730.17 


728.4 ± 0.6 


724.5 


17 


3 


775.69 


770.39 




765.6 


18 


3 


810.81 


798.35 


810.5 ± 0.4 


805.1 


19 


3 


850.78 


853.81 


849.9 ± 0.7 


846.1 


20 


3 


891.24 


894.42 




888.7 


21 


3 


930.84 


933.81 






22 


2 


962.48 


965.38 


971.7 ± 0.6 





a Wc do not list the 1 — 1.2 MOST modes until they are confirmed. 
b The observational uncertainty is ±0.40/j.Hz (quoted from original work). 
c The observational uncertainty is ±0.44/^Hz (quoted from original work). 
d Modcs used in QDG search. 

identified as radial order n + 1 in original work. 
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